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Abstract 

In this paper, we investigate a general integrable nonlocal coupled nonlinear schrodinger 
(NLS) system with the the parity-time (PT) symmetry, which contains not only the nonlocal 
self-phase modulation and the nonlocal cross-phase modulation, but also the nonlocal four-wave 
mixing terms. This nonlocal coupled NLS system is a nonlocal version of a coupled NLS system. 
The general Nth Darboux transformation for the nonlocal coupled NLS equation is constructed. 
By using the Darboux transformation, its soliton solutions are obtained. 


1 Introduction 

According to the work of Bender and Boettcher [T], the PT symmetry plays a crucial role in the 
spectrum of the Hamiltonian. They proved that a wide class of non-Herimitan Hamiltons with PT 
symmetry have real and positive spectrum. This motivated the interest for many researchers on the 
PT symmetry in quantum mechanics [2-5]. For a non-Hermitian Hamiltonian H = + V (x), if V (x) 

satisfies U(x) = then it is PT symmetric. Under this condition, the Schrodinger equation 

i^t = is PT symmetric. It has been shown that the optics can provide a good condition for 
testing the theory of PT symmetry or observing the phenomenon when the PT symmetry is breaking 
[6-9]. 

As is well known, the nonlinear Schrodinger equation (NLS), 

iqt{x, t) = Qxxix, t) ± 2\q{x, t)\‘^q{x, t), (1) 

which is PT symmetric, received a wide study since the appearance of the work by Zakharov and 
Shabat [TO]. The NLS equation has many important physical applications, such as nonlinear optics 

‘Corresponding author. Email: znzhu@sjtu.edu.cn 


1 



mi, water wave [12], and plasma physics. Very recently, Ablowitz and Musslimani investigated a 
nonlocal NLS equation m- 


iqt{x,t) = qxx{x,t) ±2q{x,t)q*{-x,t)q{x,t), (2) 

which is derived from a new symmetry reduction of the well-known AKNS system, where q{x,t) is a 
complex valued function of the real variables x and t and * denotes complex conjugation. The nonlocal 
NLS equation ([2|) is an integrable system possessing the Lax pair, infinitely many conservation laws and 
it is solvable by using the inverse scattering transformation m- Like the classical NLS equation, the 
nonlocal NLS equation ([2|) is PT symmetric and can be regarded as a mathematical model describing 
wave propagation phenomena in PT symmetric nonlinear media. The nonlocal NLS equation possesses 
some new behaviors, e.g., it simultaneously admits both the bright and the dark soliton solutions |14j . 
In addition, the dark and antidark soliton interactions have been studied by use of the Darboux 
transformation method m- 

On the other hand, integrable or nonintegrable coupled NLS equations have attracted widespread 
attention. In physical area, many problems, such as the nonlinear light propagation in a birefringent 
optical fiber [T6| , the two surface waves in deep water m, Bose-Einstein condensates [T8] , are governed 

by the coupled NLS equations. In Ref. [19], integrable properties for the following general coupled 

NLS equation are investigated, 

ipt+Pxx + 2[a|pp -I- c|gp -k bpq* + b*qp*]p = 0, 
iqt + qxx + 2[a|pp -k c|gp -k bpq* + b*qp*]q = 0, 

which contains the cross-phase modulation, the self-phase modulation, and the four-wave mixing terms. 
When a = c and 5 = 0, equation ([3]) reduces to Manakov equation m- When a = —c and 5 = 0, it 
reduces to the mixed coupled NLS equation [21]. In [19], the N-soliton solutions in this system are 
obtained by the Riemann-Hilbert method. The collision of the two soltions in this system is analyzed. 
The dark soliton, breather and rogue wave of this system are given by the Hirota method and Darboux 
transformation |22l [23l [23] . 

In Ref. [25], periodic and hyperbolic soliton solutions for a number of nonlocal nonlinear equations 
have been gained. Those nonlocal nonlinear equations include nonlocal NLS equation ([2]) and the 
nonlocal coupled NLS equation, 

ipt{x,t) +pxxix,t) -k [ap{x,t)p*{-x,t) -k bq{x,t)q*{-x,t)]p{x,t) = 0, 

iqt{x,t) + qxx{x,t) + [fp{x,t)p*{-x,t) + eq{x,t)q*{-x,t)]q{x,t) = 0, (4) 

where a, b, /, e are arbitrary real numbers. In the case o = 5 = e = /ora = / = —e = —5, this system 
reduces respectively to nonlocal Manakov equation or nonlocal Mikhailov-Zakharov-Schulman (MZS) 
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equation. The authors of the Ref. [25] asked a question: whether like the nonlocal NTS equation, 
the coupled nonlocal NTS equation of Manakov type or MZS type is also an integrable system? They 
conjectured that indeed both the nonlocal Manakov and the nonlocal MZS are integrable systems. 

In this paper, motivated by the investigation for the general coupled NTS equation dSj) and the nonlocal 
coupled NTS equation (|H), we consider the following general coupled nonlocal NTS equation: 

ipt +Pxx + 2[app*{-x,t) + cqq*{-x,t) + bpq*{-x,t) + b*qp*{-x,t)]p = 0 , 

iqt + Qxx + 2[app*{-x,t) + cqq*{-x,t) + bpq*{-x,t) + b*qp*{-x,t)]q = 0, (5) 

where p = p{x,t),q = q{x,t) and the parameters a and c are real and b is complex. It is obvious 
that equation ([5|) is a nonlocal version for the coupled NTS equation ([3|). The equation ([5|) contains 
the nonlocal cross-phase modulation, the nonlocal self-phase modulation, and the nonlocal four-wave 
mixing terms. In this paper, we will demonstrate the integrability for the coupled nonlocal NTS Q 
by proposing its Lax pair. Inspired the work of the construction of Darboux transformation for the 
coupled NLS equation in Refs. [261127] . we will construct the Darboux transformation for the nonlocal 
coupled NLS equation dSj). As an application of the Darboux transformation, we will derive N-soliton 
solutions of the nonlocal coupled NLS equation ([5|). 


2 Darboux transformation for the coupled nonlocal NLS equation 



In this section, we first demonstrate the integrability of the coupled nonlocal NLS ([S]). Then we 
construct its Darboux transformation and soliton solutions. Inspired the form of the Lax pair for the 
coupled NLS equation m, the Lax pair of equation ([5|) is given by 

= = {iXa3 + P)$, 

= -{2iX^a3 + 2XP + i{P^ + Px)(73)^, ( 6 ) 

where 



f 0 p q \ 


/ I 0 0 \ 


1 

0 

0 

, 0-3 = 

0 -I 0 


1 

to 

0 

0 


iv 0 0 -1 / 


and ri = ap*{—x,t) -|- bq*{—x,t),r 2 = b*p*{—x,t) + cq*{—x,t). In fact, note that 1 T 3 P = —Pua, and 
(TsP^iTs = P^ , we can check that the zero-curvature equation Ut — Vx + [U,V] = 0 yields 

iPt - Pxx(T 3 + 2P^a3 = 0. (7) 
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Equation (7) yields the coupled nonlocal NLS ([5]). This implies the nonlocal coupled NTS equation ([5]) 
is an integrable system. Thus, both the nonlocal Manakov equation and the nonlocal MZS equation 
are integrable systems. 

Next, using a similar method of constructing the Darboux transformation for the coupled NLS 
equation ES1E7], we construct the Darboux transformation for the coupled nonlocal NLS equation 
Suppose = <h(x, t. A), $1 = ‘h(x,t,Ai) are special vector solutions for system ([6]) at A and Ai 
respectively. We define the binary function 


D($i,$) = 




( 8 ) 


i(A + A^) 

where f is the conjugate transpose, and Re{\i) / 0. If Re{Xi) = 0, we choose <hi such that <I>i satisfies 
d>J(—x,t)M<I>i = 0 , and define D(<I>i,<I>i) by 


D(d>i, $i) = liin + e, 

X —^Al Z(^A 


(9) 


where 0 is a proper constant, and 

^ ac — | 6 p 0 0 

M = 0 —cb 

V 0 h* -a 

By a direct calculation, we have a^M = Mas, {—x,t)M = MP, Px{—x,t)M = —MP^^Pcr^ = 
—a^P. Then we can verify the following equations: 

= ^U-x,t)(iX*a 3 - PU-x,t))M = <^>U-x,t)M(iX*a 3 - P), 
ax 

j^^\-xp)M = ^\-xp) (2iX*^a3 - 2X*P^{-x,t) + ia 3 {P^'^{-x,t) + P^{-x,t))^ M 
= (-X, t)M (2iA*V3 - 2A*P + i{P^ + PjfJs) . 

Thus, D(d>i,<I>) satishes 

$) = 2(A* - A)d>J(-X, t)Mu 3 $ - (-X, t)MP^. 

Introducing the following transformations for the eigenfunctions and the potentials: 

<hiD(^>i, <I>) 


d>[l] = $- 
P[1] = P+[ct3, 


D(d>i,$i) ’ 
D($i,d>i) 


( 10 ) 


4 







then we can show 


$[!], = (aa3 + P[l])<h[l], 

mt = - (2iAV3 + 2AP[1] + i(p2[i] + P,[l])cT3) ^>[1]. (11) 


This means that the Darboux transformation for the nonlocal conpled NLS equation (5) is constructed. 


Here, we describe the main idea of the proof for equation (fTT]l . Set S 
calculation yields 




A direct 


= 2iRe{Xi)S, 


SasS 




SPS 




By using those properties of the matrix S, we can verify the validity for equation (jlljl . 

Further, following the method of the constructing N-fold Darboux transformation for the coupled 
NLS in Refs. [261127j . we can obtain the N-fold Darboux transformation for the nonlocal coupled NLS 
equation ([5]). 

Suppose (i=l,2...,N) are the N different eigenfunctions for the spectrum problem (6) at A = A* 
respectively, then the N-fold Darboux transformation can be written as 

$[iV] = $ - RW-^n, (12) 


where R = (<Li, <^ 2 ,..., and 



/ D($i,$i) 

D($i,$2) • 

■ D($i,$iv) ^ 


/ D($i,4>) \ 


IV = 

D($2,^i) 

D($2,^2) • 

Q(^2,^n) 

,D = 

D($2,^) 

(13) 


V n(Vjv,^i) 


■ j 


V D($7V,4') ) 



The relation between new potential and the old one is given by 


P[N] = P+[a3,RW-^R\-x,t)M]. 


(14) 


Here we give the proof for equation = {iXcr^ + P[A^])4>[A^]. It is similar to show that the 

eigenfunction ^[N] satisfies the time evolution equation. In fact, a direct calculation gives 

$[Ar]^ = (iAcj3 + P)$ - R^W-^Q - RW-^Q^ + RLF"^LF,1F"^D 

= (iAcj 3 + P)$ - {iasRA + PR)VF“^D - RW-^R^{-x, t)Ma 3 ^ + RW-^R^{-x, t)Ma 3 RW-^<^ 

= {iXa3 + P)<^- PRW-^Q - RW-^R^ (-x, t)Ma3<^ + RW-^R^ (-x, t)Ma3RW-^^ 

- ia3RAW-^n + a3RWR\-x, t)M<^ - ia3RW-^{XI + A*)D 

= (iAcT3 + P) 4 > - PRW-^n - RW-^R^ (-x, t)Ma3^ + RW-^R^ M(T3 + a3RWR^ (-x, t)M^ 

- iXasRW-^n - a3RW-^{R^{-x, t)MR - iWA)W-^n 
= (iAa3 + P[Ar])4.[Ar], 
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where A = diag(Ai, A 2 ,..., Xn), and we have used the two identities R^{—x, t)MR = z(TyA+A*tF),and 
i{XI + A*)0 - R\-x, t)M$ = 0. 


3 Soliton solution for the nonlocal coupled NLS 

Suppose Re{Xi) / 0, and = (/i,/ 2 ,/ 3 )^ is the eigenfunction for the spectral problem ( 6 ) corre¬ 
sponding to A = Ai, then the Darboux transformation yields 

p[l] =p[ 0 ] - (c/i/|(-x,t) - b*fif 3 {-x,t)), 

4ziie(Ai) 

g[l] = q[0] - - (a/ 1 / 3 (-x,t) - bfif 2 {-x,t)), 

0 

where g = {ac - \b\‘^)fif^{-x,t) - c/ 2 / 2 *(-x,t) - af 3 f^{-x,t) + bf 3 f^{-x,t) + b*f 2 fs{-x,t). 

Suppose p = g = p^^i{k 2 x+u] 2 t)^ where uji,ki and pi are real numbers, are the seed 

solutions of nonlocal coupled NLS ([5]), then Wj and ki must satisfy 


This last equation yields ki = k2 = 0 and wi = a;2 = w. We thus see that p = pie^‘^^,q = p2e*^* are 
the seed solutions of nonlocal coupled NLS (5), where io = 2m, and m = apf + cp| + 2bFiPiP2- In order 
to solve spectral equation ( 6 ), we make the gauge transformation 

$ = L» = diag(l, 


then the spectral equation ( 6 ) leads to 




(17) 


where 


V = 


( A -ipx -ip2 ^ 

U = i(api -I- bp2) -A 0 

\ i{b*pi + CP 2 ) 0 -X j 

( -2}? + {ap\ -b cp 2 + 2bRpip2) i2Xpi 

-i2X{api + bp 2 ) 2A^ - {apf + bpip 2 ) + ui 

\ -i2X{b*pi + CP 2 ) -{b*p\ + CP 1 P 2 ) 


i2Xp2 

-{apip2 + bpl) 

2 A^ - {cpl + h*P1P2) + w 


( 18 ) 


Consider the characteristic equation of the matrix U: 


det{Tl — U) = 0, 
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(19) 


i.e., 


(r + A)(r^ — m — A^) = 0. 

It is obvious that the characteristic equation has three roots —A, 6, —6, where 6 = —y/ A^ + m. 

Case 1. The characteristic equation has three different roots and Xr ^ 0. 

In this case, we obtain the eigenfunction = {cj)i, (j) 2 ,4>3)'^ = D'i/, where T = HNv, with v = 
(1, /I, a)'^, and 

N{x,t) = 

/ 1 1 0 \ 

i{api + hp 2 ) i{api + hp 2 ) 


H = 


\ 


5 + X -5 + \ 

i{h*pi + cp2) i{h*pi + cp2) 


-P2 


5 + X 


-5 + X 


Pi 


/ 


Further, we have 




(t>2{X,5) = i{api + hp 2 ){- 


03(A,5) = i{b* Pi +cp2){ 


P(&x-{^+2\&)t) 

A + (5 

pi[&x-(^+2\S)t) 


+ / 3 - 


+ P 


^—i{Sx+(^—2X5)t) 

A-(i ^ 

^—i{Sx+{^ —2X5)t) 




A+5 ' A-5 

and chi = (/i,/2,/3)^ = (</>i(Ai,<5i),,/.2(Ai,<5i),</.3(Ai,5i))^. 

Note that ^^-x,t)M^ = {-x,t)H^MDHN = N\-x,t)H'iMHN. Set B = ( 5 ij) 3 x 3 = 

H^MH, then we have 

/ 2Xr 2Xr ^ \ 


B = {ac — | 6 p) 


Ar + 5r Xr — i5i 
2 Xr 2Xr 


V 


Ar + i5i Xr — 5r 
0 0 


0 

0 

m 


| 6 P — ac / 


We thus obtain 


<^^-X,t)M^ = 5^^ed(5+5*)c<;-2(A5-A*5*)i) ^^^^^-i((S-S+x+2{X5+X-S*)t) ^ ^*Q^^^ia5-S+x-2{X5+X-5*)t) 


Here we focus on the case of /3 7 ^ 0, Ai E R, and Af + m < 0. We can show that in this case 
<h|(—X,t)M<hi 7 ^ 0 for arbitrary x and t. This means that nonsigular and local solution for nonlocal 
coupled NLS equation ([5]) can be obtained. We denote — m = —in, then 

= 2{ac - \h\^)2fK, 
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where 


—2k,x /Q* \ ^2kx 




Xi + in ' Xi — IK ' '^' 2(ac —|6p) 

In the following, for the sake of simplicity, we set A = Ai. By a direct calculation, we obtain 

2^^Kix-2Xt) 


p[l] = pie-‘[l - 


^^g-.(x-2At)) -.(x+2At) ^*^n(x-2M) ^ + ^ 

IH X + in X — iti pi(ac—|6p) ’ 

( 20 ) 


„r-,i „ 2A(e-(-2^*)+/3e-(-2^*)),e-(-+2^*) , /3*e-(-+2A*) , ia*(V2 + api) „_,^,+^ 2 a 2 +.H)m 

= -^ 

( 21 ) 

We remark here that if set a = 0, the solution yields the one of the nonlocal NLS equation, which has 
been given [El- 

Next we will analyse the asymptonic behavior of the solutions (|20ll and (1211) . We do this for solution 
p[l], since the asymptonic analysis for the solution g[l] is similar. 

We rewrite the solution p[l] as the following two different forms, respectively, 

’Ai X + iK X — in pi(ac —|op) 


where 


IHx = 1 + 


j^X^-2K{x-2Xt) i^*^^2K{x+2Xt) ^|^|2g-2ax+4AKt 


X + iK 


+ 


X — iK 


+ dr 6°""'' - 


2(ac- |6|2) 


+ 


P^^—2{K,x-\-2Xt) ^^2K{x—2Xt) ^ I ^,1 2^—22Ax— 4/vAt 


X + ZK 


+ 


X — IK 


+ W- 


m\a\ e 


2(ac-|6|2) 


( 22 ) 

„r-|l _ T r^tri 2 A(e 2 ^(^ 2^(x+2At) *a*(cp2 + &>l) i(Ax+(2A2 + ^)iWx+2A0u 

-^ X + ZK +A^" pi(ac-|6|2) " 

(23) 


2'! ’ 


(A) Along the line x + 2Xt = 0, we have 

_L / 3 *e^"r+ 2 At) _ ia+cp 2 +b*pi) ^-i{Xx+{2X‘^ + %)t)++x+2Xt)\' 

A A+ 2 K ' A— 2 K n-\ (ar.— \b\‘^) a 


p[l] = 1 - 


pi{ac-\b[^) 

-1 I 

' X—in 


= Pie 


Au)t 


- -J^tanhiOA + + ^Vi) i(Ax+(2AW^)t)-£i-< 

X + in A + di(ac-|6|2) e 


sech{6i) 


Xt 


—oo 


( 24 ) 



































p[l] = 1 - 


oX/Q/' e 2^(x+2At) _ ify* (cp2+fo*pi) K{x-\- 2 Xt) \ '' 

XA-ik. ' X—i.K "fTTfar^^lhi^ / 


A+2k; ' X—iK pi{ac—\b\^) 

/3Ae“2{«^a3+2Ai) 

A+ 2 K 


l/3l^ 


= /Oie 




-A 


where 


X — IK X — iK p* pi[ac — \b\^) 


Oi = k{x + 2\t) + ei + iC, 2ei = /n( ), 2C = Arg{^-^)-, 

\X — IK\ a — IK, 


\t —^ “hex) 


(25) 


02 = k{x + 2Xt) +£2 + iC, -2£2 = ln{ 


|A| 


\/3*{X + iK)\- 


(B) Along the line x — 2At = 0, we have 


p[l] = 1 - 


p[l] = 1 - 


,,, , 2A(l + /3e-2-(-2Ai)) A _ 


where 


X + iK + /3Ae 

2 A(e 2 '^(^- 2 At) + 

)^^2K(x-2\t) _j_ 


= Pie 


—A iK , / N 1 
^—I- ' , . tanh{pi) , 


X + iK X + iK 


= pie 


ibjt 


-X 


IK 


-—h - —tanh{p 2 ) , 


X — iK X — iK 


P 2 = k(x - 2Xt) + P 2 + ii', 2p2 = ln{ 


|A| 


\/3{X-iK)\- 


Xt —)• —oo 
(26) 

Xt —y Aoo 
(27) 


= k{x - 2Xt) + pi+ iv, -2pi = /n( ), -2v = Arg{ ), 

\X + ik\ X + IK 


We can see that the phase shifts along the line x + 2Xt = 0 and the line x — 2Xt = 0 are equal, 
since |ei — 62 ] = \pi — P 2 \ = It follows from (l26l) and ([27)1 that the extreme value of |p[l]| 

II |A 4 ^|m(^^ when Xsin{2iy) > Kcos{2iy), the solution p[l] is an antidark soliton. Note that 

sin{2i') = ’ eos{2v) = we can get the conclusion that the soliton along the line 

x — 2Xt = 0 is antidark when A/3/ < 0, and is dark when A/3/ > 0. Similarly, when a*{cp 2 + b*pi) = 0, 
along the line x + 2Xt = 0, the soliton is dark if XIm{/3{X — iK)"^) < 0, and is antidark if XIm{/3{X — 
iK)"^) > 0. Those phenomenons are the same as the solution for the nonlocal NLS equation. But, 
when a*{cp 2 + b*pi) 7 ^ 0, there are some differences along the line x + 2Xt = 0. In fact, it follows from 
the equations (1241) and ([25p that the solution p[l] is a combination of the bright and the dark soliton 
or the bright and the antidark soliton. Especially, when pi = 0,a*{cp2 + b*pi) / 0, the solution p[l] 
degenerates to a bright soliton with a phase shift. The solutions of the nonlocal coupled NLS equation 
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(5) are more richer than the ones of the nonlocal NLS equation. 

Case 2. iie(Ai) = 0 and ap\ + C/O^ + ‘ibpipip 2 + Af > 0. 

In this case, to obtain the solution p\l\ and (/[I], we must carefully determine <I>i). Choosing 

/ 1 \ / p*(<5a;+(f- 

i{api + bp 2 ) 


^ = D 


5 + \ 

i{h* Pi + cp2) 

\ b + \ J 


^i{5x+{^—2XS)t) _ 


( gi(<5a;+(f-2A(5)t^ \ 

i{api + hp2) nsx-{%+2\5)i) 

5 + \ 

i{b*pi + CP2) j(Sx-(^+2X5)t) 

s + x 


then we have 

^U- 


Note that, for the different spectral A and Ai, we have A^ — 5^ = Af — and 

api + cp2 + 2bRpip2 Ai + (5i 


(A + ,5)((5i -Ai) 


A + 5 


= 0 . 


We thus obtain 


+ + J((5+5^)x-(2X5+2X^5^)t) _ 

^ ’ (A + (I)((5 + di) 




— ac) 


^2iSix-\-ASikit 


Si 


0, 


where 6 = h|^L_lf£) ^ The solutions p[l\ and (^[1] are given by 


p[l] = Pie 
g[l] = p2e"‘ 


(5i + Ai 


+ 


25i 


1 


Ai - Si 5i- Ai e2*'5ia;+45ifcib + 1 


iujt f Si + Xi 


+ 


2Si 


1 


Al - Si di - Ai e2*'5ia;+45ifcii) _g 1 


(28) 


= {ac - |6|2) (1 - e^iis+s,)x-i2XS+2X,5,)t)_ (29) 

V (A + ojfOl - All / 


(30) 


This means that $t(-x,t)M$i = 0, and lim ^i( x,t)M^ ^ t{^ — ^^ 2 iSi(x- 2 Xit)_ Finally, 

A^>Ai i(A — Al) 0i 

n(<I>i,d>i) is determined by 


(31) 


We remark here that the solution is singular at {x,t) = (^^^^^,0),/c G N. 

Case 3. The characteristic equation has two multiple roots 5 = 0. 

In this case, let m < 0, then Ai = y/—m. We obtain the eigenfunction <I>i = (/i, / 2 , /s)^, where 


h 

f2 

fs 


(1 + /3(x - 2Ait))e 2 , 


ijapi + bp2) 
Al 

i{b*pi + cp2) 
Al 


(l + £+/3(a:- 
(l + £ + /3(:t 


2Ait)^ e 2 ^^-ap 2 e^^ Aix+ 2 Aft)^ 
- 2 Ait)^ 
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Thus, solutions p\l\ and (/[I] of the nonlocal coupled NLS equation are given by 


2Ait) + 1)(1 - ^ - I5*{x + 2Ait) + ^ 

Ai ^(ac-|6|^)pi J 


p[l]=pie-* (^l--(/3(x- 
<l[ll = P2e“‘ (l - |(/3(a: - 2Ait) + 1)(1 + 2Ait) - 

where = (1 — j3*{x + 2Ait))(l + /3(x - 2Ait)) + {I - — /3*(x + 2Ait))(l + ^ + /3(x - 2Ait)) - 


I 12 

g- 2 iAix of a = 0 yields the rational solution of the nonlocal NLS equation (5). 


{ac-|fe|^) 


4 2-soliton solutions 


By using N-fold Darboux transformation, we can obtain the general wave solution for nonlocal coupled 
NLS equation ([5]). Suppose = (/n,/ 2 i,/si)^ (i=l)2...,N) are the N different eigenfunctions for the 
spectrum problem (6) at A = A, corresponding to the nonzero background p[0] = pie*‘^*,g[0] = 
N-soliton solution is given by 


p[N] = - 2FiW-^{cFl{-x,t) - h*Fl{-x,t)), 

g[Ar] = - 2FiW-^{aFl{-x,t) - 6Fj(-x,t)), 


where 

-^1 = (/ii, /i2, • • • , Iin), 

F 2 = (/ 2 I, / 22 , • • • , f 2 N), 

F3 = {fsi, /32, • • • , Isn)- 

By use of the relation of = D'kj, D = diag{\, , 6“*^^*), we have the determinant representation 
of the solution 


^r/vrl ^ detjW/2—X) _ (V \ _ 

~ det{W/2) ~ K^ijjNxN — ^ 


«|JV1 = r = (r«)„xK = e 


—iu)t 


(cFj(-x,t) - b*Fl{-x,t))Fi/pi, 
[aFl{-x,t) - hFl{-x,t))Fi/p 2 . 


(33) 


Let us give the case of = 2 in detail. Set the spectral parameters Ai and A 2 are two real numbers. 
Then we have = (/ifc, f 2 k, fzkf, /c = 1, 2, where 

= g*(42:+(f-2Afe4)t) _|_ ^^g-i(4a;-(^+2Afc4)L^ 

^ i(api +fep2) ^»(4a;-(^+2A;.4)f) I + V2) ^-2(4x+(^-2A,4)f) _ 

+ Afc Afc — 5k 

^ ijb*Pi + CP 2 ) ^i(S,x-(^+ 2 XMt) I + +a^p^e-*(Afc^- 2 A 2 i)_ 

bk + Afc Afc — Sk 
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iKk- We thus have 


Case 1 A| + m < 0, /c = 1, 2. 

In this case, 6^ = = — 


^ \ •/I7|2 4AiKit p 2K.1X ^*g2Kix |2g4AiKit oiXix 

0 ($i, $i) = z(| 6|2 - ac)(-^ 

Ai Ai + ihi\ Ai — 2 Ki Ai 2ai 


0 ($i,$ 2 ) = 2 i(| 6 | 2 -ac)( 


g(K2—2(A2K2+AiKi)t l^^g—{K2+Ki)x—2{\iKi—\2K2)t 


Al + A 2 + i{Kl — K 2 ) 


+ 


Al + A 2 + i{Ki + K 2 ) 


^ ^.g(.2+M)x-2(A2K2-AiKi)i ^ ^.^^g(M-K2)a7+2(AiM+A27C2)t ^ . ma2«t +A,)x-2(A^-A?W. 

Al + A 2 — i(Ki + K 2 ) Al + A 2 — ^(^1 ~ K 2 ) Al + A 2 


0 ($ 2 ,$i) = 2 i(| 6 | 2 -ac)( 


g(Kl—K 2 )a:— 2 (A 2 K 2 +AiKi)t 


Al + A 2 + i(K2 - Ki) 


+ 


Al + A 2 + *(^^1 + K 2 ) 


^ ^.g(7C2+.Ox-2(AiKi-A2.2h ^ ;g|^^g(.2-M)a.+2(AiM+A2.2)t ^ ^ , mai«^ ^_,(a,+,,^.+2(A?-A^W. 

(Al + A2) — i{Ki + K2) Al + A2 — i{K -2 ~ ^1) Al + A2 


0($2,^2) = i{\b\^ - ac){ 


^ ^ ^ ^ •H^g- 2 iA 2 x. 

A 2 A 2 + 1^2 A 2 — iH2 A 2 2 A 2 


and 


, , -Kl(a;+2Ait) Ki{x+2\it) in-* 4 -/)* n, 

Xu = idiP - + '■ '>■)« ) 

Al + iKi Al — ^^^l pi[ac — \b\‘^) 

Kl(a;+2Ait) pKi{x+2\it) In'* -I-/)* ni'lp~h'^l*+2'^l4+^i) 

Xi 2 = i(W - ^-+ A-^- 7ai(cfe + tpi)e 

Ai+zivi Al-^^^l pi{ac-\b\^) ^ 

g-K 2 (a;+ 2 A 2 t) gK 2 (a;+ 2 A 2 i) ia* (c /12 + *) 


X 21 = i(| 6|2 - + /3ie-'‘i(^-2^i*))( 


A 2 + iK2 


X 2 — iK2 


/ii(ac- | 6 | 2 ) 


-K2{x+2\2t) K2{x+2\2t) 7’n'*('/-nr, -I-/)*niV“h-^2a;+2A|t+^i) 

X22 = ^(| 6 |^ - ac)(e"^("- 2 "^*) + ® + b pi)e 


A 2 + i/i 2 


A 2 - iK2 


Pi{ac-W) 


-K,i{x+2\p) Ki{x+2\it) ini*(hn^ rim\p~Al^+‘^>'h+%b) 

y,i = i(|6|2_ac)(e"i("-2^i*)+/3ie-"i("-2^i*))(^^—-+ ® ^ ia^[bp 2 + api)e 


Al + IKl 


Al — IKl 


+ 


p^(ac-m 


Al + iKi Al — ini p2[ac—\h\‘^) 

p-K 2 (x+ 2 \ 2 t) pK 2 {x+ 2 \ 2 t) 7 ’n'*('/)nr, -I- nni ip“d'’' 2 a;+ 2 A|t+^t) 

Ul = idiP - ocXe^l'-^'O + ^-+ A-— + ”'2(V2 + api)e 

A 2 + ^«:2 A 2 -IK 2 P2{ac-\b\^) 

-K2(x+2\2t) pK2{x+2\2t) 7’n'*('?)no-I-nni ip“d'’'2a;+2A|t+^t) 

^22 = ^(| 6|2 - ac)(e"^("- 2 "^*) + + a/nje 


A 2 + iK2 


A 2 — iK 2 


/ 92 (ac- | 6 | 2 ) 
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Case 2. m < 0, + m < 0, and 5i = —iki; 62 = 0, and A 2 = y^—m. 

In this case, we have 


n($i,$2) = ^ — ^(1 + ^2{x - 
z(Ai + A 2 ) 

(oc - | 6 p)m ^^^ g-/ci(a;+ 2 Ait) ^ ^*gKi(*+ 2 Ait) 


iA 2 (Ai + A 2 ) A 2 

_ i(-(Ai+A 2 )x+ 2 (A 2 -A 2 )t) 

i(Ai + A2) 


Ai + mi Ai — iKi 


11(^2, $ 1 ) = ^ -^(e«l(x-2Ait) ^ ^^g-^i(x-2Ait))^^ -mix + 2X2t)) 

z(Ai + A 2 ) 

(ac - | 6 | 2 )m ^ p^^-Ki(x- 2 Xit) ^ 


iX2{Xi + A 2 ) A 2 

_ «2«l"^ i(_(A,+A2)x+2(A?-AgW 

*(Ai+A 2 ) 


Ai — zKi Ai + iKi 


n{^ 2 , $2) = + / 32 (a: - 2 A 2 t))(l - / 32 *(x + 2 A 2 t) 


2 iA 2 

{ac — | 6 p)m 


^/32^ l« 2 pm 2 a 2 X 


2iAi + /52 (x - 2A2t) + ^)(1 - /32(^ + 2A2t) " " 2iA2 


X12 = i(|^P - ac)(l + P2{x - 2X2t)) 


'g-Kl(a;+2Ait) ^Ki{x+2\it) ia* (c /92 + t*) 


X 21 = i{\bf - ac)(e'"i(*- 2 Aii) ^ ^^g-Ki(x+ 2 Aii))^ 


Ai + iKi Ai — iKi 

-m 


/Ji(ac- | 6 | 2 ) 


Ai 


/3| r,\ \ ^ ia2(c/?2 + ^*/0l)e h^2a:+2A^t+“t) 


- :^{x + 2 X 2 t) + - - 

A 2 A 2 

X22 = i{\h\^ - ac)(l + / 32 (x - 2X2t) 


/9i(oc-|6|2) 

-*/3| /?2% , o^ , 1 ia^(cp2 + &Vi)e-'(^="+2^2t+fd' 

^--(x + 2A,i) + --|4|2,- 


Note that r2($i, <hi),Xu, Yu are same as the ones for the first case. Thus, two-soliton solution is 
given by the formula (|33p with N = 2. We remark here that if set = 0, then X = Y, and thus the 
two-soliton solution of nonlocal coupled NLS (5) yields the one of the nonlocal NLS equation. 


5 Conclusions and Remarks 

In this paper, we have investigated a nonlocal coupled NLS equation. We have demonstrated its 
integrability by proposing its Lax pair. We have constructed its Nth Darboux transformation. By 
using the Darboux transformation, we have derived its soliton solutions. The results obtained in this 
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paper including the Lax pair, the Darboux transformation, and the soliton solutions can yield the 
corresponding ones for nonlocal NLS equation. In the future, we will investigate the discrete version 
of the nonlocal coupled NLS equation (5). 
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